Quantum simulation of 3 He impurities and of 4 He interstitials in solid 4 He 
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We have studied the role of an atomic 3 He impurity and an interstitial 4 He atom in two- and 
three-dimensional solid 4 He using path integral Monte Carlo (PIMC) simulation. We find that when 
a substitutional 3 He impurity is introduced, the impurity becomes localized and occupies an ideal 
lattice site. When an interstitial 3 He impurity is introduced in the 4 He solid, we find that the 
impurity becomes localized at a substitutional position and, thus, promotes the extra 4 He atom to 
the interstitial space. As a consequence we find that the one-body density matrix (OBDM) and 
the superfluid fraction, for the case of a 4 He solid with an interstitial impurity, are very similar 
to those calculated for a 4 He solid with a 4 He interstitial atom. Namely, while the off-diagonal 
OBDM approaches zero exponentially with increasing particle displacement for the "pure" solid, an 
interstitial 4 He atom or a 3 He impurity appear to enhance it at long distances. Finally, the effective 
mass of the 3 He impurity quasiparticle in 2D and 3D crystalline 4 He is estimated. 

PACS numbers: 05.30.Jp, 67.80.B-, 67.80. bd, 67.80.dj 



I. INTRODUCTION 

The torsional oscillator experiments of Kim and 
Chan 1 ., where at low temperature a drop in the moment 
of inertia is observed, have motivated a number of com- 
putational studies^— of solid 4 He as well as various the- 
oretical proposals^— to explain the observation. 

There is evidence of a very strong dependence of the 
superfluid response on the 3 He impurity concentration 14 
as well as other well known facts^ about the role of im- 
purities in solid 4 Hei£~— . Proposals for the possible role 
of 3 He impurities in solid 4 He have a long history and 
date back in the late sixties^ and seventie a 20 ' 21 . In ad- 
dition, there are several experimental studies of the NMR 
relaxation of such impurities in solid helium 22 . It is, 
therefore, of great interest to study the role of impuri- 
ties in solid 4 He. Boninsegni et al£ have carried out a 
path integral Monte Carlo (PIMC) simulation of three- 
dimensional (3D) solid He using the worm algorithm 
and found that vacancies phase separate. Pollet et al»£ 
and Boninsegni et al4 have also used the above PIMC 
technique to show that grain boundaries in solid 4 He and 
screw dislocations lead to superfluidity. In addition, us- 
ing the same method Pollet et al£ have shown that the 
gap to create vacancies closes by applying a moderate 
pressure. 

In the present paper, motivated by the recent exper- 
imental and theoretical activity on the possible role of 
the 3 He impurities in solid 4 He, we study the role of 
a 3 He impurity and of an interstitial 4 He atom in two- 
dimensional (2D) and 3D solid 4 He using PIMC simu- 
lation. In addition to the motivation generated by the 
previous discussed experimental activity, this problem is 
of interest in its own right because it is not really known 
what happens locally when one injects a 3 He atom in 2D 
or 3D solid 4 He, and this can be studied by quantum 
simulation. In particular, we use the worm algorithm^ to 
simulate the 2D and 3D solid helium in the presence of 



such crystalline defects. We present results of the radial 
distribution functions and off-diagonal one-body density 
matrix (OBDM) for the following cases, (a) Pure solid 
4 He at somewhat above but near the liquid-solid melting 
density (a = 0.070 A~ 2 or a = 0.026 A" 3 ), (b) A single 
substitutional 3 He impurity in solid 4 He. (c) An inter- 
stitial 4 He atom (defect). This atom is identical to the 
other 4 He atoms and, therefore, it participates in permu- 
tation cycles, (d) An interstitial 3 He impurity in solid 
4 He. 

We find that an initial interstitial impurity quickly re- 
laxes to a regular lattice site of the 4 He solid by creating 
an interstitial 4 He atom as was proposed in Ref. I 111 Fur- 
thermore, we find that introducing such interstitial impu- 
rities in 4 He solid greatly enhances the long-distance part 
of the off-diagonal OBDM. This enhancement as well as 
the calculated superfluid response is comparable to that 
of interstitial 4 He atomic defects. It is quite possible that 
at a finite density interstitial 4 He atoms phase separate 
as do vacancies^. In such case the enhancement of the 
OBDM at long distances and of the superfluid density, 
due to a single interstitial 3 He or 4 He atom found in the 
present paper, may disappear when a finite density of 
such impurities or interstitials is introduced. However, 
interstitial atoms or impurity atoms can bind to already 
existing defects, such as dislocations or disclinations (es- 
pecially in 2D) and this tendency for phase separation 
may be avoided. In general, it is of great value to know 
what happens locally in the 2D and in the 3D crystalline 
4 He when a 3 He impurity or an interstitial 4 He atom is 
introduced. 

The paper is organized as follows. In Sec.|H]we briefly 
describe the PIMC method used to study this system. 
In Sec. IIIH we present and discuss the pair distribution 
functions for the case of a 2D and 3D solid with and 
without the introduction of a 3 He impurity and He in- 
terstitial atom. The energetics of creating such atomic 
defects in the 2D and 3D solids as well as the calcula- 
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tion of the effective mass of 3 He impurity in solid 4 He 
is discussed in Sec. |TV] In Sec. [V] we present the results 
for the one-body density matrix, the superfluid density 
and a histogram of the number of particles involved in 
the same permutation cycle for the cases (a-d) above for 
2D and 3D solid helium. Finally, the main findings as 
well as the limitations of the present work are discussed 
in Sec. EH 

II. SIMULATION DETAILS 

Using an approximation for the density matrix that is 
accurate to fourth order— in r, we use 320 time slices to 
reach a simulation temperature of IK. We have collected 
data from 2500 continuous iterations for our simulations 
in 2D, and ^1000 continuous iterations for our simula- 
tions in 3D. Each iteration consists of 500 Monte Carlo 
moves. 

All simulated atoms considered in our present stud- 
ies are isotopes of helium and therefore interact via the 
same potential. We use the Azia^i potential to model 
both the 4 He- 4 He interaction and the 3 He- 4 He interac- 
tion. With the exception of the 3 He impurity atom the 
rest are all 4 He atoms which will be treated appropriately 
to simulate their bosonic nature. The impurity atom is 
distinguishable from the "background" 4 He atoms. 

Our simulation cell is designed to accommodate cither 
a 2D 56-site triangular lattice that is very nearly square 
(25.86 A x 25.60 A), or a 3D 180-site hexagonal close- 
packed lattice (18.35 A x 19.07 Ax 17.98 A). In both 
cases of the 2D and 3D lattices we have used periodic 
boundary conditions. The density of lattice sites is fixed 
at 0.0846A- 2 (2D) and 0.0286A" 3 (3D). We will use the 
term pure solid for the case where there is exactly one 
He atom per lattice site. The term substitutional solid 
will be used when a single 4 He atom is removed from 
the pure solid and is replaced with an impurity atom. 
Additionally, the term interstitial solid will be used when 
a single atom (either 4 He or an impurity) is added to the 
pure solid. 

III. DISTRIBUTION FUNCTIONS 

A. Two-dimensional solid 4 He 

How does the impurity atom affect the pair distribu- 
tion function (744 of the 4 He atoms of the underlying 
solid? We find that when a substitutional impurity is 
introduced it becomes localized and occupies an ideal lat- 
tice position with its own zero-point motion determined 
by its different mass. In Fig. Q] we present the calcu- 
lated (744 (r) radial distribution for pairs of 4 He atoms for 
the four different case systems studied: (a) pure solid 
4 He (dashed line) (b) the 4 He solid with a substitutional 
3 He impurity (also dashed line) (c) the 4 He solid with an 
interstitial 4 He defect (solid line), and (d) the 4 He solid 
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FIG. 1: Here we show the radial distribution function, (744(7-), 
for pairs of 4 He atoms in two dimensions. The organizational 
structure of the 4 He atoms does not change in the presence 
of a substitutional impurity. However, when an interstitial 
defect or impurity is present, we can see that 344(f) becomes 
less peaked at the nearest-neighbor distance lattice positions. 
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FIG. 2: A snapshot of a space-time configuration for the 2D 
triangular solid, after thermalization, and starting from a con- 
figuration with an interstitial 3 He atom. Each atom's trajec- 
tory in imaginary time appears as fractal covering a finite 
size spot. The crosses (red in the online version) are the 4 He 
atoms and the circles (blue in the online version) is the 3 He 
impurity atom. 



with an interstitial 3 He impurity (also solid line). Within 
the accuracy of our results we cannot discern any differ- 
ence in the (744 distribution function for the cases of the 
pure solid and the substitutional impurities. When an 
interstitial impurity is present in the 4 He solid, we find 
that the impurity becomes localized at a substitutional 
position, thereby promoting the extra 4 He atom to the 
interstitial band. This is shown by the snapshot space- 
time configuration shown in Fig. [2 Notice that while the 
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FIG. 3: Top: The pair distribution function 344(2:, y) for pure 
2D solid 4 He. Middle: Contour plot of the distribution func- 
tion for pairs of 4 He atoms, gn(x,y) in the pure solid. Bot- 
tom: the difference Sga between 344 of the pure (solid) and 
that of the solid with interstitial impurity. 



initial configuration has an interstitial 3 He impurity, in 
the configuration obtained after thermalization (shown in 
Fig. [5]) the 3 He becomes substitutional by promoting an 
interstitial 4 He atom. Namely, in the equilibrium config- 
uration, shown in Fig.[5J the 3 He atom, in our lattice with 
periodic boundary conditions, is located in a regular tri- 
angular lattice position surrounded by six 4 He atoms. In 




FIG. 4: Top panel: Contour plot 544(2;, y) for the interstitial 
solid. This function is independent of the type of defect or 
impurity. Bottom panel: Contour plot of the distribution 
function 5134(2;, y) for the same interstitial solid. 



addition, a 4 He atom has been promoted to the intersti- 
tial space which creates larger density fluctuations in the 
crystalline arrangement in some parts of the system. As 
a consequence of this fact ga(r), in Fig. [3l is less peaked 
at the lattice positions. In Fig. |3j (top) the calculated 
pair distribution function 1744(2;,?/) for pure 2D solid He 
is shown and in Fig. [3J (middle) we present the contour 
plot of the same (744 (x,y). This function is nearly iden- 
tical for the substitutional solid (which is not shown, as 
it looks exactly alike). This implies that the introduced 
substitutional impurity becomes localized and it only af- 
fects its neighboring atoms. In the case of an interstitial 
impurity the difference in the 344 distribution function, 
as discussed above and shown in Fig. Q] and Fig. [3J (bot- 
tom), is significant because the added impurity takes the 
position of a 4 He atom and, thus, there is an extra 4 He 
atom that necessarily becomes interstitial. In the bottom 
panel of Fig. [3] we plot 5g^(x,y), the difference between 
#44(2;, y) of the pure solid and the solid with a single in- 
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FIG. 5: The radial distribution function for pairs of He 
atoms in the three-dimensional HCP lattice simulation cell, 
344 (r). The organizational structure of the 4 He atoms does 
not change in the presence of a substitutional impurity. How- 
ever, when an interstitial defect or impurity is present, by 
looking at Ag4i(r) (scale on the right) we can see that ga(r) 
becomes less peaked at the nearest-neighbor distance lattice 
positions. 



really does have a 4 He atom in the interstitial space, 
and also that the interstitial 3 He solid has relaxed into 
a space where the 3 He interstitial atom has become sub- 
stitutional, and in doing so promoted a 4 He atom to the 
interstitial band. 
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terstitial impurity. Notice that the extra atom is truly 
interstitial since the 544 is reduced by an amount in the 
neighborhood of the ideal lattice positions and enhanced 
in the interstitial space by the same amount. It was ver- 
ified through integration in the enhanced regions (or the 
reduced regions) finding exactly one extra 4 He atom in 
the interstitial regions. 

Our finding that the interstitial impurity becomes lo- 
calized at regular lattice sites can be further illustrated by 
comparing the contour plots of the 344(2, y) and 1734(2;, y) 
for the case where we have a 4 He solid with an intersti- 
tial impurity. In the top panel of Fig. [4j we present the 
contour plot of the distribution function 1744(2;, y) for the 
case of a 4 He solid with an interstitial impurity. Within 
the accuracy of the discretization of the probability den- 
sity of the contour plot this function is independent of 
the type of defect or impurity. In the lower panel is the 
distribution function 534(2;, y) for pairs consisting of the 
impurity atom and one 4 He atom. Because the contour 
plots for both 334(2;, y) and 344(2;, y) are identical in shape 
and in form, we may surmise that the impurity atoms are 
located at lattice sites. 

B. Three-dimensional solid 4 He 

In Fig. [5] we show 344(r) for the 3D system. As in 2D, 
we find that the pure solid and the substitutional solid 
are nearly identical in structure, as are the two inter- 
stitial solids. Also shown is the difference, Ag^r), be- 
tween 344 (r) the pure solid and the interstitial solid. As 
expected, 344 (r) for both interstitial solids is less peaked 
at lattice positions compared to the pure and substitu- 
tional solids. This indicates that the 4 He interstitial solid 



FIG. 6: Potential energy of a 3 He atom placed either substitu- 
tionally solid line (blue in the online version) or interstitially 
dashed line (green in the online version) into triangular solid 
4 He. After a brief relaxation, both energy values remain close 
except for occasional "blips" in the potential energy of the 
initially interstitial 3 He atom. 




FIG. 7: Snapshot of a space-time configuration of the He 
interstitial solid in a "blip" of elevated potential energy that 
appears after thermalization. Red crosses represent the 4 He 
atoms at each imaginary time slice, while blue circles repre- 
sent the 3 He impurity. The 3 He atom can be seen to be at a 
region of local disorder. 



5 



IV. ENERGETICS OF IMPURITY AND 
INTERSTITIAL 

If a 3 He atom, initially placed in the interstitial re- 
gion of a triangular solid of 4 He atoms, relaxes onto a 
lattice site by promotion of a 4 He atom to the intersti- 
tial space, this should be seen in the energy values of the 
simulated atoms. In Fig. [6] we show the potential energy 
of a 3 He atom in the substitutional and interstitial 3 He 
solids. A short relaxation time can be seen for the in- 
terstitial solid, as the 3 He atom relaxes onto the lattice. 
After that, the potential energy of a 3 He atom in both 
systems is almost the same. After 600 iterations, a small 
bump is seen in the energy of the (initially) interstitial 
3 He atom. A snapshot of the atomic configuration at this 
elevated energy value is shown in Fig. [7] The 3 He atom 
is no longer at an equilibrium lattice position, but rather 
at what appears to be a possible edge dislocation. This 
is not entirely unexpected, as a 3 He atom in solid 4 He 
exhibits a high rate of diffusion. Such "blips" in the en- 
ergy of the 3 He in the interstitial solid occur occasionally 
throughout our simulation, but account for no more than 
5% of configurations. 



Int. 
Int. 



2D 3D 
4 He - Pure 4 He 50.27 ± 0.54 K 22.4 ± 1.3 K 
3 He - Sub. 3 He 50.41 ± 0.55 K 24.1 ± 1.2 K 



TABLE I: Excitation energy of an interstitial 4 He atom, as 
calculated by the difference in energy between (1) the pure 
solid and the interstitial 4 He solid, and (2) the substitutional 
solid and the interstitial 3 He solid. 

In Table U we show the activation energy of an intersti- 
tial 4 He atomic defect. This is calculated by subtracting 
the total energy of the pure solid from the total energy of 
the interstitial 4 He solid. If the interstitial 3 He solid is ac- 
tually the substitutional solid with an added interstitial 
4 He atom, as we propose it is based on the distribution 
functions above, then the activation energy can also be 
calculated by subtracting the total energy of the substi- 
tutional solid from the total energy of the interstitial 3 He 
solid. We find that both methods give activation energies 
in agreement with one another. 

We have also estimated the effective mass of the 3 He 
impurity in solid 2D and 3D 4 He using our data on 
the imaginary time diffusion following Ref. [25|. Namely, 
we approximate the low-energy (which dominates the 
long time evolution) impurity quasiparticle spectrum by 
the dispersion near the T point of the Brillouin zone of 
both the triangular 2D solid and of the hexagonal closed 
packed 3D lattice 



E(k) 



h 2 k 2 

2m* ' 



(1) 



It is straightforward to carry out the imaginary-time evo- 
lution for this spectrum and to calculate the average of 
(r(0) — r(r)) 2 , where r(r) is the impurity coordinate in 
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FIG. 8: Top: The ratio m/m* as a function of r for 2D and 
T = 1 K (Solid line) and T = 0.5K (dashed line). Bottom: 
The ratio m/m* as a function of r for 3D and T = 1 K for 
substitutional and interstitial 3 He impurities. 



imaginary time. We find that, 



m 

= hm 

m* t->/3/2 



((r(O)-r(r)) 2 ) 



2d\ 



-iP-r) 



(2) 



where A = Ti 2 / (2m) and d is the dimensionality. In Fig. [5] 
we plot the right-hand-side of the above equation as cal- 
culated from our simulation for the 2D (Fig. Ettop)) and 
3D (Fig. ISJbottom)) case. We find that in the 2D case 
the effective mass ratio of the 3 He impurity at T = IK is 
5.10 ± 0.02 while at T = 0.5 K it increases to 9.06 ± 0.04. 
In the 3D case we have available results only for T = 1 
K, where the substitutional and the interstitial impurity 
masses are found to be 5.67 ± 0.03 and 5.47 ± 0.04 re- 
spectively. 
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FIG. 9: The one-body density matrix, n(r). Although no FIG. 11: Histogram of relative frequency of accepted particle 
difference is observed between the pure solid and the substi- permutations for various number of particles in 2D. 
tutional solid, the interstitial solid clearly shows a significant 
enhancement of n(r) quantity. 
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V. OFF-DIAGONAL ONE-BODY DENSITY 
MATRIX 




FIG. 10: The one-body density matrix, n(r), of the 180-site 
HCP system in three dimensions. Although no difference is 
observed between the pure solid and the substitutional solid, 
the interstitial solid clearly shows a significant enhancement 
of n(r) quantity. 

In Fig. [9] we compare the one-body density matrix n(r) 
for (a) defect- free solid 4 He (solid line), (b) solid 4 He 
with a substitutional 3 He impurity (dotted line), (c) solid 
4 He with an interstitial 4 He defect (long-dashed line), (d) 
solid 4 He with an interstitial 3 He impurity (dashed line), 
and Notice that the substitutional 3 He impurity and the 
pure solid have similar one-body density matrices. On 
the contrary, a 4 He solid with interstitial 3 He impurity 



FIG. 12: Histogram of relative frequency of accepted particle 
permutations for various number of particles in 3D. 



and a 4 He solid with interstitial 4 He atoms have one-body 
density matrices which are significantly enhanced at long 
distances. This result agrees with the fact that winding 
numbers (and hence superflow) are observed in the inter- 
stitial solid (see Table [TTJ) . Notice that these superfluid 
fractions are very high considering that the simulation 
was carried out at IK. The reason for these high su- 

2D 3D 
Interstitial 3 He 0.021(7) 0.007(4) 
Interstitial 4 He 0.011(6) 0.012(5) 

TABLE II: Supersolid fraction, p s /p, in the presence of an 
interstitial atom. No global permutations were observed for 
the perfect lattice and the substitutional impurity. 
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perfluid fractions is finite size effects. These results for 
the superfluid fraction arc presented in order to make the 
case that a interstitial impurity has a very similar effect 
on the superfluid fraction and OBDM as an interstitial 
4 He atom. 

In Fig. [10] we compare the one-body density matrix 
for the 3D results. As in 2D, both the pure solid and 
the substitutional solid show exponential decay of n(r). 
Although the enhancement of n(r) at large distance is not 
obvious for the interstitial 3 He solid, it is very obvious for 
the interstitial 4 He solid. This may be due to a shorter 
MC run as compared to the 2D data. In any case, once 
again both interstitial solids display superfluidity, while 
the pure and substitutional solids do not (see Table III]) . 

In Fig. [UJ (for 2D) and in Fig. [TJ (for 3D) we present 
a histogram of cycles (i.e., how often in the simulation 
we encounter cycles of exchanges involving a given num- 
ber of particles). Notice that in both 2D and 3D case, 
the pure solid and the 3 He substitutional solid has only 
one or two particle permutation cycles, while when an 
interstitial 3 He or A He atom is introduced, it gives rise 
to permutations involving up to a 10 atom chain, which 
is as long as the longest possible distance in our lattice. 
This indicates that the result may not be a finite-size 
effect. 



VI. DISCUSSION 

One of the main conclusions of the present paper is 
that the added interstitial impurity in both 2D and 3D 
4 He becomes substitutional by creating a interstitial 4 He 
defect; we believe that this result is firm and it is not 
subject to finite size effects. Furthermore, we find that 
the effective mass of a 3 He impurity atom in both 2D 
and 3D solid 4 He is large at T = 1 K (m*/m ~ 5) and 
at a lower temperature of 500 mK in 2D it becomes even 
larger (m*/m ~ 9). 

In addition, we find that the above mentioned effect 
(i.e., the promotion of a 4 He atom to the interstitial band 
by the impurity) gives rise to a non-zero superfluid re- 
sponse and a significant enhancement of the OBDM at 
long-distances. This suggests that, provided that this ef- 
fect persists when a finite density of 3 He impurities is 
present and, provided that such a metastable state can 
be created and maintained, 4 He solid with such impu- 
rities should be a supersolid. However, this can not be 
established by the present calculation done for a single 
impurity in a pure 4 He solid and it depends on a num- 
ber of other factors. For example, while we have clearly 
demonstrated that a single 3 He impurity acts as a donor 
of He atoms to the interstitial ("conduction") band, the 
fate of these freed bosonic "carriers" is not certain when 
there is a finite density of 3 He impurities. In this case 
the created interstitial He atoms can phase-separate in 



a similar way as vacancies do 6 -, or they may bind to exist- 
ing defects, such as, dislocations, domain walls, or grain 
boundaries or even remain free. It is not clear that such 
interstitial defects exist in the 4 He solid caused by 3 He 
or other impurities. This is an issue which could depend 
on the process of the crystal growth^—. 

A 2D 4 He solid only exists as films on substrates, such 
as on graphite. The phase diagrams of first, second, third 
and fourth layer of 4 He on graphite, as a function cover- 
age, has been studied by PIMC simulation 26 *. The role of 
substrate corrugations, which is missing from the present 
simulation of the ideal 2D 4 He, is important and the in- 
terplay of these substrate potential corrugations with the 
helium-helium interaction gives rise to a wealth of inter- 
esting phases 26 . It is quite possible, however, that the 
main conclusion of the present paper, that introducing 
an interstitial 3 He impurity is solid 4 He leads to the pro- 
motion of a 4 He atom to an interstitial position while the 
3 He impurity becomes substitutional, may remain valid 
even in the case of substrate corrugations. 

The superfluid response which was calculated at 1 K 
and is given in Table lHl is very large considering the fact 
that the calculation was done at such a high temperature. 
This is a finite-size effect but at a much lower temper- 
ature the superfluid response is expected to be greater. 
A calculation of the superfluid density at a significantly 
lower temperature requires much larger computational 
time scales in order to be able to accurately sample it. 
In the 3D case, the zero temperature condensate frac- 
tion obtained as the asymptotic value (infinite distance 
value) of the off-diagonal OBDM at zero temperature, 
is much smaller by at least two orders of magnitude (as 
seen from Fig fTU)) . Therefore, as is well-known, there is a 
large factor relating the superfluid response and the ac- 
tual condensate fraction. It is clear that introducing just 
a single impurity and taking the infinite volume limit (or 
infinite area limit in 2D), the superfluid density and the 
condensate fraction should vanish. It is interesting, how- 
ever, the fact that the ratio of the values of both these 
two quantities to the impurity fraction (the impurity frac- 
tion is 1/N, where N is the total number of 4 He atoms 
considered) is a number of order unity. These reported 
results on the off-diagonal OBDM and superfluid density, 
have only a qualitative value and one cannot draw firm 
conclusions because of a) finite size effects and b) they 
refer to the case of a single 3 He impurity or single 4 He 
interstitial. 
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